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Abstract 


An  analytical  and  computational  study  of  the  gasification  and  oxidation  of  an  energetic 
liquid  fuel  droplet  is  presented.  Single-step,  finite-rate,  Arrhenius  reaction  rate  expressions 
are  used  for  exothermic  liquid-phase  decomposition  and  gas-phase  oxidation.  The  liquid  fuel 
is  assumed  to  decompose  to  a  gaseous  product  at  a  fixed  number  of  bubble  sites  per  unit 
mass  (specified  a  priori )  within  the  droplet.  Decomposed  gas  escapes  the  droplet  surface 
by:  (1)  decomposition  (gasification)  of  the  droplet  surface,  (2)  decomposition  at  the  surface 
of  bubbles  that  connect  with  the  droplet  surface,  and  (3)  escape  of  gas  inside  bubbles  due 
to  droplet  surface  regression.  Without  oxidation,  results  are  compared  between  one  model 
wherein  gaseous  fuel  leaves  the  droplet  due  to  decomposition  of  the  droplet  surface  and 
bubbles  that  connect  with  the  droplet  surface  when  the  void  fraction  exceeds  a  critical  value 
{4>c)~  and  another  model  wherein  the  droplet  mass  decreases  due  to  discontinuous  bubble 
bursting  at  the  droplet  surface.  The  models  agree  well  (with  the  exception  of  oscillations  in 
the  droplet  radius  predicted  by  the  latter  model)  in  the  limit  of  <pc  1. 

The  transient,  two-phase,  governing  equations  are  solved  numerically  for  various  values  of 
the  nondimensional  reaction  rate  coefficients  (for  both  decomposition  and  oxidation),  heats 
of  decomposition  and  oxidation,  number  of  bubbles  per  unit  mass  ( N/m ),  and  ambient 
temperature  and  pressure.  Consistent  with  simplified  scaling  for  the  limit  of  chemical  rate 
control,  the  droplet  lifetime  i]d.*  is  strongly  dependent  on  the  nondimensional  decomposition 
rate  constant  and  activation  energy,  and  less  strongly  on  the  number  of  bubbles  per  unit  mass, 
ambient  pressure,  and  heat  of  decomposition.  Increasing  the  ratio  of  gas-phase  to  liquid- 
phase  thermal  conductivities  increases  slightly.  The  droplet  lifetime  is  a  monotouically 
increasing  function  of  initial  droplet  radius  but  too  weak  for  diffusion  rate  control:  the  results 
are  closer  to  the  chemical  rate  control  limit. 

After  an  ignition  delay  period,  the  flame  radius  is  predicted  to  increase  nearly  linearly 
with  time  until  the  droplet  is  gasified.  After  this  time  (%.,),  the  flame  radius  decreases 
with  time.  The  variation  of  flame  radius  with  time  differs  from  classical  droplet  burning  due 
to  the  exothermic  decomposition  process  that  determines  the  gasification  rate.  Simplified 
scaling  previously  derived  for  the  droplet  lifetime  also  correlates  the  effect  of  decomposition 
parameters  on  the  flame  behavior. 

Gas-phase  oxidation  does  not  appreciably  affect  the  droplet  lifetime,  for  the  selected  base 
case  values  of  the  above  parameters,  because  droplet  heating  is  controlled  primarily  by  the 
liquid-phase  decomposition.  As  the  decomposition  rate  is  reduced  (e.g..  by  reducing  N/m). 
the  time  scale  for  heat  conduction  from  the  flame  to  the  droplet  becomes  comparable'  to 
that  of  liquid  decomposition,  and  hence  gas-phase  oxidation  significantly  reduces  r)(U.  For 
the  base  case,  liquid-phase  decomposition  increases  the  flame  temperature  by  approximately 
6%.  Therefore,  energetic  liquid  fuels  can  be  expected  to  vaporize  and  burn  more  rapidly  and 
to  yield  more  thermal  energy  than  conventional  hydrocarbon  fuels. 


NOMENCLATURE 

Aig  surface  area  of  liquid/gas  interface  per  unit  mass  of  fluid 
Cp  specific  heat 


1 


Ec  activation  energy  for  gas-phase  oxidation 

Ed  activation  energy  for  liquid-phase  decomposition 

G  gasification  rate  per  unit  surface  area  of  liquid/gas  interface 

kc  rate  constant  for  gas-phase  oxidation 

kd  rate  constant  for  liquid-phase  decomposition 

Leg  Lewis  number,  A '/(pg,*CpV,) 

p  ambient  pressure 

Qc  heat  release  rate  per  unit  mass  of  fuel  due  to  gas-phase  oxidation 

Qd  heat  release  rate  per  unit  mass  of  fuel  due  to  liquid-phase  decomposition 

Rb  gas  pocket  (bubble)  equivalent  radius 

r  nondimensional  radius,  r' j r' 

r'0  initial  droplet  radius 

Hu  universal  gas  constant 

t  time 

T  temperature 

vr  radial  velocity 

Yi  mass  fraction  of  species  i 

Greek  Symbols 

P  weighting  factor  for  the  reduction  of  bubble  surface  area  due  to  bubble  connec  ting 
77  nondimensional  time,  7/XL/rf 

77 Cj*  time  at  which  99%  of  the  initial  droplet  mass  has  oxidized  in  the  gas  phase* 

77d,*  time  at  which  99%  of  the  initial  droplet  mass  has  gasified  (droplet  lifetime) 

ujcd  species  consumption  rate  due  to  gas-phase  oxidation 
u>d  gasification  (decomposition)  rate  per  unit  volume  of  fluid 
4>  void  (gas  volume)  fraction 

<f>c  critical  packing  value;  bubbles  connect  when  p>  4>c 

mass  coordinate  divided  by  47 r 
gas-phase  reference  density  (air  at  P'x,  T'x) 
p  material  gas  density 

Pi  material  liquid  density 

p  bulk  gas  density 

Pt  bulk  liquid  density 

r  time  in  von  Mises  plane 

6  nondimensional  temperature.  ( T '  —  —  T'0),  where  T'0  =  300  K 

£  normalized  mass  coordinate.  i)/p's 

Subscripts 

*  reference  property  evaluated  at  300  K  for  the  liquid  phase  and  at  for 
the  gas  phase 

fl  flame  (point  of  maximum  fuel  consumption  rate) 
g  gas  phase;  subscript  sometimes  omitted 

i  gas-phase  species  (e.g.,  F,  fuel  vapor;  P.  products) 

/  liquid  phase 

s  droplet  surface 
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Other  nomenclature  is  defined  when  first  used  in  the  text.  Dimensional  variables  are  denoted 
with  a  prime  (/)  and  nondimensional  variables  are  unprimed. 

INTRODUCTION 

The  additional  flight  range  available  with  high-energy,  strained-molecule  fuels  such  as 
benzvalene  is  highly  attractive.  The  relative  chemical  instability  associated  with  the  new 
fuels  introduces  new  physics  and  chemistry  into  the  droplet  vaporization  problem.  The 
exploration  of  this  new  science  is  important  and  challenging  from  both  the  practical  and  the 
intellectual  points  of  view.  The  most  critical  new  feature  is  the  occurrence  of  exothermic 
reactions  in  the  liquid  phase.  Upon  heating,  some  of  the  chemical  bonding  in  the  liquid  fuel 
will  be  broken,  releasing  more  energy  and  yielding  simpler  hydrocarbons.  Thermodynamic 
analysis  indicates  that  the  decomposition  energy  per  unit  mass  of  a  strained  hydrocarbon 
fuel  is  one  order  of  magnitude  less  than  the  energy  per  unit  mass  of  fuel  released  in  oxidation. 
In  oxidation,  however,  the  energy  is  shared  with  a  mass  one  order  of  magnitude  larger  than 
the  mass  of  fuel  due  to  mixing  of  the  fuel  with  air.  Adiabatic  temperatures  for  oxidation 
of  conventional  fuels  and  for  decomposition  of  strained  fuels  are  therefore  of  the  same  order 
of  magnitude.  For  example,  according  to  Moriarty  and  Rao  (1993).  benzvalene  has  a  strain 
energy  of  4.13  x  106  J/kg.  Published  heats  of  combustion  and  vaporization  for  benzene  are 
4.06  x  107  J/kg  and  3.95  x  105  J/kg,  respectively. 

Law  and  co-workers  (Lee  et  al..  1988,  1989  and  1992)  conducted  experimental  studies 
of  the  vaporization,  combustion,  microexplosion,  and  thermophysical  and  thermochemical 
properties  of  organic  azides.  They  found  that  the  quiescent  gasification  rate  of  droplets  of 
organic  azides  is  greater  than  the  corresponding  hydrocarbons,  and  they  attributed  this  to 
liquid-phase  decomposition  either  within  the  droplet  or  at  the  droplet  surface  (Lee  et  ai. 
1988).  While  monoazides  were  found  to  closely  fit  a  Clausius-Clapeyron  relation  over  a  broad 
temperature  range,  diazides  heavier  than  diazidohexane  did  not  follow  such  a  relation.  It 
was  hypothesized  that  decomposition  both  in  the  liquid  phase  as  well  as  upon  gasification 
can  occur  in  the  heavy  diazides.  This  hypothesis  was  later  supported  by  droplet  experi¬ 
ments  performed  at  elevated  pressures  (Lee  et  al.,  1992).  As  the  pressure  was  increased 
in  these  experiments,  a  large  gasification  rate  increase  was  observed  for  cliazides  but  not 
for  monoazides  or  hydrocarbons  of  equal  volatility.  Liquid-phase  reactions  were  found  to 
increase  due  to  dihalide  addition  and  droplet  temperature  elevation.  By  using  the  burning 
rate  constant  in  the  classical  d2-Law  expression  to  interpret  their  experimental  results,  Lee 
et  al.  (1989)  conjectured  that  decomposition  occurs  in  the  droplet  interior  and/or  at  tin1 
droplet  surface  rather  than  in  the  surrounding  gas  phase.  Lee  et  al.  (1989).  however,  did 
not  observe  evidence  of  liquid  decomposition  in  monoazide  droplets  heated  up  to  500  K. 
Microexplosions  were  observed  in  organic  azide  droplets  (Lee  et  al..  1988),  which  is  a  strong 
indicator  of  bubble  formation.  It  is  possible  that  (1)  the  bubbles  were  too  small  to  detect, 
and/or  (2)  the  microexplosion  of  the  droplets  was  caused  by  bubbles  formed  so  rapidly  that 
they  were  not  detected. 

It  is  interesting  to  note  that  liquid-phase  reactions  might  also  occur  in  the  combus¬ 
tion  of  high-energy  solid  fuels  such  as  HMX  and  boron/poly(BAMO/XMMO)  fuel-rich  solid 
propellants.  Palopoli  and  Brill  (1991)  describe  subsurface  reactions  that  create  a  heteroge¬ 
neous  “foam  zone”  in  the  condensed  phase  of  HMX.  A  subsurface  reaction  zone  was  also 
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described  by  Hsieh  et  al.  (1991)  for  the  combustion  of  boron/poly(BAMO/NMMO).  The 
liquid  layer  that  is  observed  during  the  rapid  thermal  decomposition  of  HMX  is  formed 
by  condensed-phase  liquefaction  with  concomitant  decomposition  rather  than  true  melting 
(Brill  and  Karpowicz,  1982;  Palopoli  and  Brill.  1991). 

Experiments  have  not  yet  revealed  useful  information  on  the  rates  of  fuel  decomposition, 
liquid  vaporization,  and  gas-phase  mixing.  In  particular,  droplet  experiments  have  only  been 
performed  with  the  high-energy  fuels  strongly  diluted  bv  conventional  fuels  (Law.  1995). 
Theoretical  and  numerical  analyses  should  reveal  vital  information  about  the  reaction  zone 
in  the  liquid  where  reactants  and  gaseous  products  exist  together. 

The  high  energy  fuel  can  conceivably  decompose  in  at  least  three  ways:  (1)  high  tem¬ 
peratures  might  cause  the  strained  fuel  to  vaporize  and  then  decompose  and  burn  in  the 
gas  phase;  (2)  the  liquid  decomposes  to  a  liquid  product  that  gasifies:  and  (3)  the  liquid 
decomposes  to  a  gaseous  product.  The  first  path  requires  a  greater  stability  for  the  strained 
fuel  in  order  for  the  fuel  to  vaporize  without  decomposing.  The  heavier  molecule  of  the 
original  strained  fuel  would  be  short-lived  under  gas-phase  collisions.  The  very  large  heat 
release  due  to  decomposition  would  probably  make  the  lifetime  of  the  product  as  a  liquid 
negligible,  thus  making  path  (2)  unlikely.  This  leaves  the  third  path  as  a  likely  case  and  the 
one  that  is  studied  herein. 

The  third  path  requires  a  very  hot  liquid  vaporizing  in  depth.  On  the  fine  scale,  gasifi¬ 
cation  has  a  distribution  of  nucleation  sites  throughout  the  hot  liquid  interior.  At  the  high 
temperatures,  nucleation  is  occurring  faster  than  outer  droplet  surface  regression  so  that 
bubbles  should  result.  The  liquid  is  being  heated  so  rapidly  that  the  outer  surface  cannot 
accommodate  vaporization  quickly  enough  to  maintain  phase  equilibrium.  As  the  liquid 
superheats,  other  nucleation  sites  in  the  liquid  interior  are  given  time  to  develop. 

The  liquid-phase  reaction  is  most  likely  to  occur  in  the  warmer  portion  of  the  droplet  near 
the  surface  rather  than  the  cooler  domain  in  the  center  of  the  droplet.  The  exothermic  liquid- 
phase  reaction  will  augment  the  heating  of  the  reactants  (strained  fuel)  and  the  products 
(still  a  hydrocarbon  fuel)  in  this  reacting  domain.  Existing  droplet  vaporization  analyses  do 
not  treat  the  change  of  phase  within  the  droplet  interior  that  can  result  from  this  heat  release. 
It  is  not  reasonable  to  expect  the  liquid-phase  reaction  to  occur  in  a  very  thin  spherical  shell 
around  the  edge  of  the  droplet;  this  would  require  very  fast  reactions  at  the  relatively  low 
liquid  temperatures.  More  likely,  the  reaction  occurs  over  a  finite  domain  where  both  liquid 
and  gas  exist.  With  large  enough  gaseous  content,  all  of  the  gas  cannot  remain  in  solution, 
but  rather  bubbles  should  form.  It  is  reasonable  to  expect,  therefore,  that  a  bubbly  region 
will  exist  near  the  surface  that  consists  of  liquid  fuel  reactant  and  gaseous  fuel  products. 
This  two-phase  bubbly  region  exists  between  the  pure  liquid  and  the  surrounding  gas. 

Figures  1  and  2  show  schematics  of  the  vaporizing  liquid-fuel  energetic  droplet  with  and 
without  gas-phase  oxidation.  The  model  both  temporally  and  spatially  resolves  all  three 
regions  of  the  problem:  the  surrounding  gas  film,  two-phase  bubbly  layer,  and  the  liquid 
droplet  core.  The  decomposed  energetic  liquid  fuel  forms  bubbles  in  the  droplet,  causing  tin* 
droplet  to  swell  before  the  droplet  radius  decreases  due  to  the  transport  of  decomposed  liquid 
to  the  surrounding  gas  film.  Fuel  vapor  is  assumed  to  cross  the  interface  of  the  droplet  and 
enter  the  gas  phase  (i.e.,  gasify)  by  three  mechanisms:  (1)  decomposition  (gasification)  of 
the  droplet  surface,  (2)  decomposition  at  the  surface  of  bubbles  in  the  droplet  that  connect 
with  the  droplet  surface  either  by  intersecting  the  droplet  surface  or  by  connecting  with 
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Figure  1:  Schematic  of  the  vaporizing  liquid-fuel  energetic  droplet  without  gas-phase  oxida¬ 
tion:  (a)  shows  the  initial  droplet  (t i  =  0).  while  (b)  and  (c)  show  the  droplet  at  progressively 
larger  times,  to  and  £3. 


/  Two-phase 


Figure  2:  Schematic  of  the  vaporizing  liquid- fuel  energetic  droplet  with  gas-phase  oxidation. 
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other  bubbles  that  intersect  the  droplet  surface,  and  (3)  escape  of  gas  inside  a  bubble  due 
to  droplet  surface  regression.  Results  using  this  decomposition  model  are  compared  below 
to  an  alternative  model  that  assumes  that  the  droplet  mass  decreases  due  to  discontinuous 
bubble  bursting  at  the  droplet  surface.  Unless  otherwise  specified,  bubbles  are  assumed  to 
connect  with  the  droplet  surface  (mechanism  #2  above)  when  the  void  fraction  exceeds  a 
critical  packing  value,  (pc  =  0.8. 

The  model’s  input  parameters  include  initial  droplet  temperature,  ambient  pressure  and 
temperature,  gas  and  liquid  densities,  specific  heat,  Lewis  number,  and  the  following  nondi- 
mensional  parameters:  rate  constants  (for  both  the  decomposition  rate  and  oxidation  rate 
expressions),  heats  of  decomposition  and  oxidation,  number  of  bubbles  per  unit  mass,  and 
ratio  of  gas-phase  to  liquid-phase  thermal  conductivities.  A  parametric  study  is  performed 
below  on  some  of  these  key  parameters  to  determine  the  effects  of  gas-phase  oxidation  on  the 
droplet  lifetime,  as  well  as  the  effects  of  the  liquid-phase  decomposition  on  the  flame  radius, 
flame  temperature,  and  gas-phase  profiles  of  temperature  and  species  mass  fractions. 


ANALYSIS 

The  continuity,  energy,  and  species  conservation  equations  are  solved  numerically  using 
a  von  Mises  formulation  with  a  mass  coordinate  that  is  normalized  by  the  droplet  mass. 
Gas-phase  oxidation  and  liquid-phase  decomposition  are  modelled  by  exothermic,  single- 
step,  finite-rate  Arrhenius  chemical  reaction  rate  expressions.  In  the  simplified  model  used 
here,  mass  diffusion  in  the  liquid  phase  is  not  considered.  Therefore,  the  gas  bubbles  in  the 
liquid  droplet  contain  only  decomposed  fuel  vapor  that  does  not  dissolve  back  into  the  liquid 
phase.  The  following  additional  assumptions  are  made:  (1)  mass  diffusion  in  the  liquid  phase, 
transport  of  heat  by  mass  diffusion,  kinetic  energy,  and  viscous  dissipation  are  neglected: 
(2)  the  momentum  equation  is  replaced  by  assuming  zero  pressure  gradients;  (3)  spherically- 
symmetric  geometry,  (4)  bubbles  are  spherical  and  do  not  migrate;  (5)  the  specific  heats  of 
the  vapor  and  the  liquid  are  assumed  equal  and  constant  (a  reasonable  assumption  for  higher 
molecular  weight  hydrocarbons);  (6)  the  decomposed  gas  is  assumed  to  be  benzene  (il//.-  = 
78.11);  (7)  droplet  shattering,  microexplosions  and  the  breakup  of  bubbles  are  neglected 
(the  number  of  bubbles  per  unit  mass  is  assumed  constant):  (8)  the  change  in  gas/liquid 
interface  area  due  to  the  coalescence  of  bubbles  is  neglected:  and  (9)  the  material  density 
of  the  liquid  is  constant,  p\  =  1000  kg/m3.  This  value  of  liquid  density  is  representative 
of  many  energetic  liquid  fuels  that  commonly  exist  in  liquid  form.  Varying  p\  from  700  to 
1300  kg/m3  does  not  change  qualitatively  the  results  of  this  study.  For  comparison.  Law  and 
co-workers  (Law,  1996)  have  studied  liquid  dihydrobenzvalene  ( p\  =  854  kg/m3)  and  liquid 
methylcubane  ( p\  =  1010  kg/m3). 

The  key  assumptions  of  the  model  are:  (1)  the  use  of  single-step  reaction  rate  expressions 
for  decomposition  and  gas-phase  oxidation,  and  (2)  the  a  priori  specification  of  bubble  site 
locations.  Our  objective  is  to  study  qualitatively  the  effects  of  various  fundamental  nondi- 
mensional  groupings  on  the  gasification  and  oxidation  of  an  energetic  liquid  fuel.  The  use 
of  one-step  decomposition  kinetics  clearly  oversimplifies  the  complicated  chemistry  involved 
in  the  problem  but  is  both  practical  and  necessary  without  additional  information  about 
the  individual  decomposition  reaction  steps.  Gas-phase  oxidation  is  expected  to  occur  rel¬ 
atively  fast,  so  the  use  of  a  one-step  finite-rate  reaction  expression  is  arguably  a  sufficient 
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improvement  to  an  assumption  of  infinite-rate  reactions.  The  present  model  assumes  there 
are  a  specified  number  of  bubble  sites  per  unit  mass  (N/m)  rather  than  relying  on  kinetic 
theory  to  predict  the  location  of  bubble  sites.  In  practice,  it  is  possible  that  the  number  <»f 
bubbles  might  be  controlled  through  seeding  with  tiny  particles  and  the  encouragement  <>f 
heterogeneous  nucleation. 

Let  p  denote  the  void  mass  fraction  in  a  two-phase  gas-liquid  region.  Then.  0  <  p  <  1: 
<t>—  1  in  the  gas  phase;  and  0= 0  in  the  liquid  phase.  The  bulk  density  for  the  gas  is  ji  =  op' . 
and  for  the  liquid  it  is  p't  =  (1  — 0)p{.  Under  the  assumptions  above,  the  transient,  dimensional 
conservation  equations  are  as  follows: 

Conservation  of  mass: 


d{p'  +  p'i)  ,  1  d 


dt' 


+  r'2  dr' 


,2 


(fi  +  P'l)vr]  ~  0 


(1) 


Conservation  of  energy: 

[dV  ,dT' 

(P  +  Pl^  (  W  +  Vr 


dr' 


P'  ^  Q'd  ,  Qc  _i _ 

*c'FC'p  C'pr a  dr' 


C'pdt'  +  WdC'p 


r’2\' 


2,0 


dr 

p  dr' 


C2) 


Conservation  of  species  mass  in  the  gas  phase  surrounding  the  droplet  (where  0  =  1’ 


'd\\  ,  ,d\ 


i  d 


j  l^*i,  r^*i\  x  ^  /  .2  isr,  v  1 1  l  / 

P  (  Oi<  +  r  Or  I  )  “  ,./2  Or!  1  P  Or!  : 


dYi 


dt' 


dr 


r'2  dr' 


dr 1 


Here  u'c  i  is  the  consumption  rate  of  species  i  due  to  gas-phase  oxidation. 

The  property  variations  (in  space  and  time)  result  from  varying  temperature  and  com¬ 
position  of  the  two-phase  zone.  Several  studies  are  available  to  predict  thermal  conductivity 
of  a  mixture.  Here,  we  use  the  general  mixture  rule  from  Nielsen  (1978)  with  the  reduced 
concentration  term  equal  to  unity: 


_  1  +  ABo 
^  1  -  Bo 

where  A  is  a  constant  that  depends  on  the  shape  of  the  bubbles,  state  of  agglomeration,  and 
nature  of  the  interface.  We  assume  spherical  bubbles  in  our  case  so  that  .4  =  2  (Nielsen. 
1978).  The  constant  B  is  given  by 


o  _  m  - 1 
A'/ A;  +  .4 

Note  that  =  \\  when  0  =  0  and  \'2  <p  =  X'g  when  o  =  1. 
The  equation  of  state  is  given  by 

P  =  4>p'  =  4>p'M/KT ' 


{■)) 


(G) 


where  M  is  the  mean  molecular  weight  of  the  mixture  of  air  and  decomposed  gas. 


The  gasification  rate  (or  liquid  mass  loss  rate)  per  unit  area  of  gas-liquid  interface  is 
assumed  to  be 

G' =  k'^M-E'JKT')  [kg'/in2-  s]  (7) 

The  rate  constant  k'd  and  the  activation  energy  E'd  depend  upon  the  type  of  the  liquid  fuel. 
The  problem  is  closed  for  a  particular  value  of  k'(i  and  E'd.  If  A'ltJ  denotes  the  total  surface 
area  of  liquid/gas  interface  per  unit  mass  of  fluid,  then  the  total  gasification  rati'  per  unit 
volume  of  fluid  is 

u'd  =  +  p'i)G'  [kg/m3  •  s]  (8) 

If  R'b  denotes  the  equivalent  average  bubble  radius,  then  A'lg  =  Xr:R'b{N'  j  in')  in  the  interior 
of  the  droplet  if  the  bubbles  do  not  touch  each  other.  As  discussed  below,  different  equations 
for  A\g  arise  due  to  either  gasification  of  the  droplet  surface  or  the  reduction  of  bubble  surface 
area  due  to  bubble  connecting. 

The  number  of  bubbles  (bubble  centers)  per  unit  mass,  N'/w1.  is  related  to  the  number 
of  bubbles  per  unit  volume  (the  bubble  number  density,  n')  by 


The  void  fraction  is 


.V 

m' 


— - -  =  constant 

p'  +  p'i 


,  47 rgv 
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Also,  the  rate  of  change  of  bubble  radius  is  given  by 

dr1  p'g  Zpfg  dr' 


(9) 

(10) 


(11) 


where  p'g  is  the  density  of  the  gaseous  products  from  the  decomposition  of  the  liquid.  Eq.  (11) 
is  derived  from  a  Lagrangian  analysis  (see  von  Mises  transformation  below). 

The  fuel  consumption  rate  due  to  gas-phase  oxidation  is 


=  Kp^’fYo,  exp {-E'JU'J')  (12) 

Assuming  a  stoichiometric  reaction  of  benzene,  'Jc  0.2  =  3.077u;'  F  and  a/  P  =  —4.077a;'. 
For  the  present  study,  values  of  k'c  =  3.6  x  109  m3/kg-s  and  E[  =  30  kcal/mole  were  si'lected 
from  hydrocarbon  fuel  data  (Westbrook  and  Dryer.  1981),  and  Q'c  =  4.06  x  10'  J/kg  is  that 
of  benzene. 


Nondimensional  Equations  After  Coordinate  Transformations 

When  the  droplet  swells  due  to  bubble  formation,  the  droplet  radius  changes  whereas  the 
droplet  mass  remains  fixed  if  gasification  of  the  droplet  surface  were  negligible.  Therefore,  a 
von  Mises  transformation  (Chervinskv,  1969;  Bhatia  and  Sirignano.  1992)  is  applied  so  that 
coordinates  are  changed  from  an  Eulerian  time  and  radius  ( t r')  to  a  Lagrangian  time  and 
mass  coordinate  (r',t/;'),  where  r'  =  t', 


gd_ 

dr'  dt'  1  r  dr' 


(13) 
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and  the  stream  function  ip'  is  defined  as 


dip' 

dt' 


-r'\p' +  p\)v'r 


dl''  -  *2,-/.- ^ 

Tr 7  -  <P  P,) 


d-n 


Note  that  Eq.  (14)  satisfies  Eq.  (1). 

In  order  to  provide  a  fixed  boundary  with  respect  to  the  regressing  droplet,  a  final 
coordinate  transformation  is  made  from  ( r',ip ')  to  (//'.£)  where 


£  =  fp'/^'s  ,  v'  =  r' , 


d 

dr' 


Vj  ^ 


d> f 


(1-3) 


Here  ip's  and  v's  are  the  droplet  mass  and  time  derivative  of  the  droplet  mass  (both  divided 
by  47 r),  respectively.  The  evaluation  of  ip's  is  discussed  below.  Note  that  ^  is  nondimensional. 
The  radius  of  any  given  point  is  determined  by: 


Wo  p’  +  p', 


(10) 


where  £  is  a  dummy  variable. 

The  nondimensional  variables  are  defined  as  follows: 


n  =  n'r'g,  m 


m 

V../.3  ’ 

POo 


T  f '  r'  P  r' 2  I'D' ' 


i.i 


v  — 


3d 
Pi7  o 


pi  _  p'i 

P  =  7-  Pi  =  =  1  -  Ps 

Pi  Pi 


p'i'1''  VPP  P  dCUTU-Tl) 


Tr  —  T7 

0  =  -— W  W 


A' 


T1/  _  T "7 

J  DO  1  O 


kr  = 


WP* 

k'r' 


V'o 


,  Qc  = 


pW 

Q'c 


<5rf  = 


Qrf 


-  T'0) 


p\v:  ^  C'p(T^  -T0) 

Other  nondimensional  groupings  important  to  this  problem  include: 


T'jT'o,  N/m,  A'/A;. 


E' 


H'T' ’  11' T' 


K 


( 1 7  > 


Using  the  above  coordinate  transformations  and  nondimensional  variables,  the  mixture 
continuity,  species  and  energy  equations  may  be  written  as 


3  r2 

lps 


C P  +  Pi? 


dvr 

af 


"(P  +  P<K 
r 


d{p  +  P;) 

dn 


fl  Wp'> 


d}]  AVYf  _  9pg,.  d  f  4  u;c,t 

dp  W  d£  ~  ip]  dt  \  9p9  dt,  )  pg 


(IS! 


(19i 


9 


dp 


Thde_ 

**  i’s  d£ 


9/V» 


Leg 


d 

dd 


AlgQdgex  p 


p  dd 
(p  +  Pi) 


v$  d(j> 
vs 


—r2(p  + Pi)  vr 
ips 


do 


(20) 


Eq.  (18)  is  solved  only  within  the  droplet  to  obtain  values  of  vr  for  the  pressure  source  term 
in  Eq.  (20).  Eq.  (19)  is  solved  only  for  the  gas  film  surrounding  the  droplet  because  liquid- 
phase  mass  diffusion  is  neglected.  The  evaluation  of  Aig  is  discussed  below.  0  =  1  outside 
of  the  droplet  because  no  bubbles  exist  there.  Hence  the  equation  for  0  in  the  gas  film 
surrounding  the  droplet  does  not  have  the  last  two  source  terms  in  Eq.  (20).  Note  also  that 
the  product  pg,*Leg  appears  rather  than  the  term  Le  used  by  Schiller  et  al.  (1996)  because 
p'gV  is  assumed  constant. 

The  nondimensional  fuel  consumption  rate  due  to  gas-phase  oxidation  is 


u c,f  =  kcp2YFY0,  exp {-E'jn'uT')  (21) 

where  kc  —  k'cp\r'2 /Vt.  The  base  case  values  of  the  nondimensional  oxidation  coefficients 
are  kc  =  2.16  x  1010  (r'Q  —  1  mm  and  £>,  evaluated  at  1000  K),  E[i'R,'uT'0  =  50.4.  and 
Qc  =  Q'c/[C'p(T^  -  T')]  =  58. 

The  rate  of  change  of  bubble  radius  is  given  by 

dRb  G  Rb  dpg  ij>a  ( dRb  Rb  dpg\ 
dp  pg  3  pg  dp  V  di  3  pg  d£  ) 


where 


G  =  kd  exp  <  - 


Tl'T' 


(T  N 

1  +  9  1.0 


In  nondimensional  form,  the  number  of  bubbles  per  unit  mass,  assumed  uniform  and 
constant,  is 


N 

m 


XpK3 


constant 


m'  p  +  Pi 

Combining  equations  (9)  and  (10)  with  p  =  dp'/p'i  and  p,  =  1  -  <fi  yields 


,  /47r\  ( N\  _o._  _ . 

<’=  (t)  {m)RbiP  +  P,)=  . 


47T 

T 


m 


1  + 


/47T\  (N\  3 

(yj  ImW-d 


(23) 


The  boundary  conditions  are  dO/dE,  =  0  at  £  =  0,  and  6=1,  YF  =  dYp/dr  =  0.  and 
Yo2  —  Po2, oo  (0.23,  unless  otherwise  specified)  at  £  — »  oo. 

Composite  properties  are  used  at  the  grid  interfaces  (Patankar,  1980).  so  only  boundary 
conditions  for  1)  are  required  at  the  interface  between  the  liquid  droplet  and  the  gas  phase. 
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The  following  mass  flux  balances  are  used  for  the  interface  boundary  conditions  for  the 
species  mass  fractions: 

9i>\ _ P-W  #.  (24a, 

dr  ) s  3 {p  +  Pi),>i  dr] 

and,  for  oxygen  and  products, 


/  d¥j\  Yj 

\dr)s  3(p  +  p/)srj  dr] 


Decomposition  Model 

In  this  section,  equations  are  derived  for  Atg  and  vs  due  to  (1)  gasification  of  the  droplet 
surface,  (2)  connecting  of  bubbles,  and  (3)  escape  of  gas  inside  a  bubble  due  to  droplet  surface 
regression.  Let  (.7  denote  the  droplet  mass  loss  rate  (divided  by  4tt)  due  to  decomposition 
and  rjjg  denote  the  droplet  mass  loss  rate  from  gas  that  escapes  due  to  regression  of  the 
droplet  surface.  If  the  ratio  of  the  two  corresponding  volumetric  loss  rates  is  the  same  as  the 
ratio  of  the  gas  and  liquid  volumes,  then  0/  and  rf'g  are  related  by 


t/’o 


=  i'r. 


4>  p\ 


<Pp'i 


g  _  P_ 
‘Pi 


(23) 


and  rps  =  ipt  +  vg.  ipt  is  related  to  the  individual  liquid/gas  interface  areas  that  contribute 
to  droplet  mass  loss  by  mechanisms  (1)  and  (2)  above  by 


(20) 


where  G,  is  the  temperature-dependent  gasification  rate  for  each  region  i.  The  factor  of  3 
appears  because  the  reference  values  for  m  and  rp  are  different. 

Consider  a  spherical  layer  of  thickness  5  =  Rb  just  below  the  droplet  surface.  Bubbles  that 
lie  partially  in  this  layer  intersect  the  droplet  surface  and  gasify  directly  to  the  surrounding 
gas  film.  If  the  fraction  of  droplet  surface  area  intersected  by  the  bubbles  is  0  (evaluated  just 
below  the  droplet  surface),  then  the  droplet  surface  area  that  gasifies  is  _T t  =  4tt/-;(1  —  o). 
If  we  further  assume  that  any  polar  angle  of  intersection  on  the  bubble  is  equally  probable, 
then  on  average  an  intersecting  bubble  has  a  surface  area  of  2t rRj  that  decomposes  below 
the  droplet  surface.  The  total  decomposing  surface  area  of  all  the  bubbles  that  intersect  the 
droplet  surface  can  be  shown  to  be  A2  =  127tr^0.  If  bubbles  do  not  connect  below  this  thin 
subsurface  layer,  then  rp[  is  determined  from  Eq.  (26)  using  areas  .4t  and  .4-2  with  G  for  these 
regions  calculated  with  6S. 

Bubble  connecting  occurs  when  the  void  fraction  exceeds  a  critical  packing  value,  o,. 
For  simply  packed  spheres,  <pc  =  7t/6,  while  4>c  =  0.7405  for  hexagonal  packing.  If  the 
bubbles  are  not  spherical  in  reality,  then  larger  values  of  oc  are  possible.  If  we  assume 
that  the  total  bubble  surface  area  is  reduced  by  a  weighting  factor  3  <  1  due  to  bubble 
connecting,  then  Alg  =  4irRl(N/m)3  is  used  in  Eq.  (20)  for  grid  points  with  q  >  o,. 
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whereas  Aig  =  4irRl(N/m)  for  grid  points  with  0  <  0C.  One  exception  to  this  is  the  Hist 
subsurface  grid  below  the  droplet  surface  (i=IDROP),  where  the  number  of  bubbles  that  do 
not  intersect  the  droplet  surface  is 

For  t=IDROP,  the  chemical  source  term  in  Eq.  (20)  includes  the  heat  release  due  to  the 
decomposition  of  areas  .4i  and  .42  as  well  as  from  the  bubbles  that  do  not  intersect  the 
droplet  surface. 

If  the  bubbles  connect  with  each  other  and  with  the  droplet  surface,  then  gasification 
within  the  bubble  can  contribute  to  vt.  In  the  present  model,  once  6  reaches  oc,  all  gas  that 
decomposes  in  the  bubble  is  assumed  to  cross  the  interface  of  the  droplet  and  go  into  the  gas 
phase.  Therefore,  the  gasification  source  term  is  omitted  from  Eq.  (22)  for  all  grid  points 
with  0  >  0C,  and  the  contribution  to  0/  is  determined  from  Eq.  (26)  using  Aigj,  G ,  and 
for  each  of  these  grid  points. 

In  the  model  without  droplet  surface  decomposition,  the  governing  equations  were  solved 
in  a  (r,  0)  coordinate  system,  and  the  mass  of  the  droplet  was  assumed  to  decrease  by 
discontinuous  bubble  bursting  at  the  droplet  surface  when  0  >  1  -  e.  where  e  —  1()-7. 
(Results  were  insensitive  to  e  providing  e  was  less  than  approximately  10-5  or  10~r>.)  The 
small  number  e  was  used  to  convert  a  region  of  nearly  all  bubbles  into  gas  film.  The  grid 
points  remained  fixed  in  the  0  coordinate  system,  and  the  location  of  the  droplet  surface 
was  defined  to  be  the  outermost  point  at  which  1  —  0  >  e.  This  model  will  be  referred 
to  herein  as  the  “(r,  0)”  model,  and  the  model  with  droplet  surface  decomposition  will  be 
referred  to  as  the  “(77, 0”  model.  The  important  differences,  however,  are  in  the  physics  of 
the  decomposition  model  and  not  in  the  coordinate  transformations. 

SOLUTION  PROCEDURE 

We  assume  that  a  pure  liquid  droplet  of  initial  temperature  90  —  0  is  injected  into  air  at 
pressure  p'  and  at  temperature  9 x.  The  calculation  proceeds  in  an  iterative  manner  for  the 
dependent  variables  V),  9 ,  vr  (velocity),  and  Rb  (bubble  radius).  Other  quantities  such  as 
void  fraction,  gasification  rate,  and  densities  are  determined  as  a  function  of  these  dependent 
variables. 

The  coupled  governing  equations  are  solved  numerically  using  a  semi-implicit,  finite- 
difference  scheme  that  is  accurate  to  second  order  in  space  and  first  order  in  time  (Patankar. 
1980).  The  solution  procedure  for  each  time  step  follows.  The  gasification  rate  crossing 
the  droplet  surface,  0s,  is  calculated  and  vs  is  updated.  Eqs.  (18)— (23)  are  then  solved  to 
determine,  in  sequential  order,  Rb  and  0  at  each  point  in  the  droplet  (after  which  thermal 
conductivities  and  densities  are  updated  and  source  terms  for  the  energy  equation  are  cal¬ 
culated),  9,  (in  the  gas  film  surrounding  the  droplet),  and  vr  in  the  droplet.  The  above 
steps  are  repeated  until  all  dependent  variables  converge.  The  radius  at  each  point  is  then 
updated  using  the  nondimensional  form  of  Eq.  (16). 

The  energy  and  species  equations  are  considered  converged  if  the  greatest  change  in  the 
value  of  a  dependent  variable  from  one  iteration  to  the  next  is  less  than  0.1%  of  tin1  value  of 
the  variable.  Excellent  convergence  is  typically  obtained  in  only  two  iterations.  Numerical 
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results  differ  by  less  than  1%  when  the  time  step  and/or  mesh  size  are  varied  by  a  factor  <>f 
two. 

The  radius  of  the  outer  boundary  in  the  gas  phase  decreases  with  time  if  £(IMAX)  is 
fixed.  For  simulations  with  the  effects  of  gas-phase  oxidation,  to  avoid  the  numerical  problem 
of  having  the  flame  reach  the  outer  boundary  of  the  gas  phase,  after  each  time  step  the  ga— 
phase  grid  is  adjusted  so  that  0(IMAX)=5O  is  fixed.  This  value  is  sufficient  to  keep  the  flame 
away  from  the  outer  boundary  throughout  a  simulation.  Linear  interpolation  is  used  to  map 
the  field  variables  in  the  gas  phase  from  the  old  to  the  new  grid.  There  is  no  noticeable 
difference  in  results  between  the  model  with  this  grid  shift  and  results  without  this  grid  shift 
technique  for  the  cases  of  (1)  no  gas-phase  oxidation,  and  (2)  with  gas-phase  oxidation  for 
times  at  which  the  flame  radius  does  not  reach  the  outer  boundary. 

Simulations  which  did  not  include  oxidation  were  terminated  when  ?/  =  >jd^  (the  droplet 
lifetime,  when  999c  of  the  initial  droplet  mass  has  gasified).  For  simulations  with  oxidation, 
decomposition  was  artificially  terminated  at  rj  =  rather  than  converting  the  remaining 
liquid  to  gas  because  the  mass  coordinate  is  normalized  by  xps.  These  cases  were  terminated 
when  either  77  =  (the  time  at  which  99%  of  the  initial  droplet  mass  has  oxidized)  or  the 
total  nondimensional  oxidation  rate  (integrated  over  the  gas  phase)  decreased  to  below  10-' 
after  the  droplet  fully  decomposed.  The  latter  condition  is  met  if  the  flame  extinguishes 
before  all  the  fuel  has  oxidized. 

In  the  resolution  of  the  large  gradients  near  the  droplet  surface,  the  nonlinear  relation 
between  the  coordinate  f  and  the  radial  coordinate  mandates  the  use  of  a  nonuniform  grid 
for  f.  A  cubic-ally  increasing  mesh  £nj,  =  n|A£,  where  nk  is  the  grid  point  number,  is 
employed.  One  exception  to  this  formula  is  the  first  grid  point  outside  the  droplet,  where 
A£(IDROP+l)  =  A£(IDROP)p5(0  =  0)  is  used  to  resolve  the  initial  heat  fluxes  at  both 
sides  of  the  droplet  surface.  Input  parameters  include  the  number  of  mesh  points  in  the 
liquid  droplet.  IDROP  (which  determines  A£),  and  the  total  number  of  mesh  points.  To  get 
at  least  one  bubble  site  in  the  first  control  volume  surrounding  the  center  of  the  droplet,  one 
must  have  IDROP  >  1  +  ( N/m )C3. 

Most  of  the  simulations  were  run  with  101  grid  points  in  the  droplet.  201  grid  points 
in  the  gas  phase,  and  a  time  step  A77  =  10-3.  Results  are  insensitive  to  near  factor  of  two 
changes  in  time  step  or  number  of  grid  points.  All  runs  were  performed  on  a  DEC  Alpha 
5000/200  workstation.  A  grid  dependency  check  was  made  by  fixing  the  radius  of  the  outer 
numerical  boundary,  R(IMAX),  and  using  a  constant  value  of  A rg.  Results  (e.g..  Rn  and 
08  versus  1 7)  do  not  change  appreciably  if:  (1)  A rg  or  At;  is  decreased  by  a  factor  of  20.  or 
(2)  the  number  of  grid  points  in  the  droplet  is  decreased  by  30%. 


RESULTS  AND  DISCUSSION 

Results  of  parametric  studies  with  and  without  gas-phase  oxidation  are  presented  sepa¬ 
rately  below.  I  nless  otherwise  specified,  all  of  the  calculations  reported  herein  used  T'a  -- 
300  K,  =  1000  K,  C'p  =  103  J/(kg-K),  p'  =  1  atm,  Leg  =  0.87,  and  3=1. 

A.  Results  Without  Gas-Phase  Oxidation 

The  important  parameters  that  we  have  studied  are  kd,  Qd.  X/m.  X'  jX\.  E'd/{Tl'uT'0). 
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p',  <j)c,  and  j3.  The  effects  of  these  parameters  on  the  decomposition  rate  and  the  liquid 
temperature  distribution  are  discussed  below. 

Simplified  scaling  can  be  used  to  compare  with  the  detailed  computational  results  and 
predict  the  effect  of  various  parameters  on  the  characteristic  time  //,/..  for  the  droplet  to 
decompose.  Taking  the  chemical  rate  control  limit  in  which  increases  only  due  to  gasifi¬ 
cation,  Eq.  (22)  gives 


Rb 


Pg  1 kd exp  - 


& d 


(27) 


where  the  unknown  constant  0*  >  1  and  accounts  for  gasification  at  a  temperature  higher 
than  Tg. 

Substitution  of  Eq.  (27)  into  Eq.  (23)  and  scaling  77  -»  when  </>  — ►  1  yields 


(28) 


Eq.  (28)  can  also  be  obtained  by  scaling  dd/di]  with  the  heat  release  due  to  chemical  de¬ 
composition  and  6  —¥  Qc  as  77  — »  r]d^.  In  the  following,  the  droplet  lifetime  is  defined  as 
the  time  at  which  95%  of  the  initial  droplet  mass  has  decomposed.  In  terms  of  dimensional 
variables,  the  above  scaling  indicates  that,  in  the  chemical  limit,  ?fdt  is  independent  of  the 
initial  droplet  radius,  r'0. 

At  77  =  0,  the  droplet  loses  mass  only  due  to  gasification  of  the  droplet  surface1  at  the 
initial  rate  of  x hi  =  3kdexp[—Ed/(7l[Tg)].  The  time  rp  at  which  the  rate  of  droplet  mass  loss 
due  to  gasification  of  the  droplet  surface  equals  the  rate  of  mass  loss  due  to  decomposition  of 
bubbles  that  intersect  the  droplet  surface  can  be  obtained  by  equating  the  areas  .4[  and  ,-U 
defined  above  (see  section  entitled  "Decomposition  Model")  and  by  assuming  that  G  is  the 
same  for  these  areas.  Thus  771  is  equal  to  the  time  at  which  <p  =  1/4,  which  from  Eqs.  (23) 
and  (27)  is 

Vl  ~  _47r(_Y/m) 

If  0*  =  1  in  Eq.  (28),  then  T]i/rjd^  =  (3 pg)~1^3  ~  0.07.  For  all  the  parameters  investigated 
below,  Eq.  (29)  correctly  predicts  77!  within  a  factor  of  two  or  less.  Using  Eq.  (28)  with 
0«  =  1  tends  to  overpredict  77^,  by  about  a  factor  of  two  for  (pc  —  0.99  and  underpredict  //,/., 
by  about  a  factor  of  two  for  <t>c  =  7r/6  for  p'  =  1  atm. 


I/O 


Pg{ at  0  =  0) 
kdexp[-E'd/(TZ'uT^)} 


(29) 


Results  of  Base  Case 

The  following  parameters  were  selected  for  the  base  case:  kd  =  10-4,  Qd  -  3.  N/m  =  I()'\ 
X'/X'i  =  0.2,  p'  =  1  atm  (p  =  1.1  x  lO"4),  and  E'J{K'UT'0 )  =  4.6.  The  base  case  value  of  (),,  is 
approximately  one  order  of  magnitude  less  than  the  heat  of  combustion  for  a  hydrocarbon 
fuel.  This  value  is  representative  of  a  strained  hydrocarbon  such  as  benzvalene.  Tin1  heat 
released  due  to  liquid-phase  decomposition  in  this  case  causes  the  flame  temperature  of  an 
energetic  liquid  fuel  to  be  approximately  10%  greater  that  that  of  a  hydrocarbon.  In  tin1 
absence  of  any  definitive  data  in  the  literature,  we  selected  E'd/(R!aT'0)  —  4.6  so  that  the 
value  of  E'd  is  one  order  of  magnitude  lower  than  the  activation  energy  for  gas-plmse  kinetics 
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of  hydrocarbon  fuels  and  is  of  the  same  order  of  magnitude  as  the  product  of  M r  and  a 
typical  hydrocarbon  latent  heat  of  vaporization.  Note  that,  for  conventional  hydrocarbon?*, 
this  product  is  effectively  an  activation  energy  for  vaporization.  The  number  of  bubbles  per 
mass,  N/m,  was  selected  so  as  to  have  at  least  one  bubble  site  in  the  smallest  computational 
control  mass.  Eq.  (28)  with  the  base  case  values  and  /q  evaluated  at  G50  I\  gives  //,/.,  =  35. 
This  value  of  is  one  order  of  magnitude  less  than  the  nondimensional  droplet  lifetime 
predicted  by  the  classical  d2-Law  for  hydrocarbon  droplets  vaporizing  in  a  non-oxidizing. 
1000  K,  1  atm  environment  (Kanury,  1975). 


Thin  lines:  (t,y)  model 
Thicker  lines:  (r|,^)  model,  <|>c=0.99 


Time,  x  or  rj 

Figure  3:  Variation  of  radius  and  mass  coordinate  at  the  droplet  surface  for  the  (r,  </’)  model 
without  droplet  surface  decomposition  and  the  (rj,  £)  model  with  droplet  surface  decompo¬ 
sition  using  data  for  the  base  case:  kd  =  10~4,  Qd  =  3.  N/m  —  106,  A'/AJ  =  0.2,  //  =  1  atm. 
and  E'd/ {VJUT'0)  =  4.6.  Two  values  of  <j)c  are  shown  for  the  (/;,£)  model. 

Figure  3  shows  a  comparison  of  rs  and  tJ)s  versus  time  for  the  two  models  using  the  base 
case  parameters.  After  <f>  — *  1  with  the  (r,  ij;)  model  without  droplet  surface  decomposition, 
the  intermittent  bursting  of  bubbles  into  the  gas  film  causes  small  oscillations  of  the  droplet 
radius  as  it  intermit  tenth-  swells  (due  to  sensible  heating  of  the  bubbles)  and  shrinks  (due  to 
the  bubble  bursting).  Conversely  with  the  (77.  £)  model  with  droplet  surface  decomposition, 
the  droplet  radius  continuously  increases  with  time  until  surface  decomposition  becomes 
appreciable,  after  which  time  the  droplet  mass  and  droplet  radius  decrease  continuously 
with  time. 

Figure  4  shows  radial  profiles  of  temperature  at  different  times  obtained  from  the  (//.  £) 
model  for  the  base  case  with  (j)c  =  0.99.  The  results  are  similar  to  those  from  the  (r,  </’)  model. 
The  nondimensional  temperature  does  not  exceed  unity  until  the  void  fraction  just  below 
the  liquid  surface  exceeds  approximately  0.9.  The  energy  released  due  to  the  decomposition 
causes  the  temperature  to  peak  near  the  droplet  surface  with  #max  — >  Qd-  A  comparison  of 
the  individual  terms  in  the  energy  equation  indicates  that  the  rate  of  temperature  increase. 
dO/drj,  just  below  the  droplet  surface  is  due  primarily  to  conduction  at  //  =  10.  At  ij  =  15. 
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Figure  4:  Radial  profiles  of  temperature  at  different  times  obtained  from  the  (//.£)  model 
for  the  base  case  with  0C  =  0.99. 

conduction  and  the  heat  release  due  to  decomposition  contribute  more  or  less  equally  to  the 
rate  of  temperature  increase  just  below  the  droplet  surface.  For  i]  >  20.  the  heat  due  to 
decomposition  is  much  greater  than  the  heat  conduction  throughout  the  droplet.  The  source 
term  due  to  the  moving  boundary  is  typically  negative  for  f  <  1,  negligible  for  //  <  10.  and 
of  the  same  order  as  the  chemical  decomposition  term  in  the  region  0.5  <  <  0.9  for  //  >  20. 

At  all  times  and  radial  positions,  the  pressure  term  in  the  energy  equation  is  negligible. 


Effects  of  <pc.  3  and  p' 

As  the  critical  packing  value,  oc.  is  decreased  in  the  (//,£)  model,  the  droplet  radius 
increases  more  gradually,  the  maximum  value  of  rs  decreases,  and  the  time  for  the  droplet 
to  decompose  fully  increases.  For  example,  =  31  for  <f>c  =  0.99,  %.*  =  49  for  0,.  =  0.8. 
and  rjd^  —  85  for  0C  =  7r/6.  These  effects  are  due  to  the  reduction  in  the  rate  of  increase  of 
Rb  and  hence  0  when  0>  0C  (see  Fig.  5)  because  the  mass  of  gas  in  the  bubbles  is  assumed 
fixed  (gas  that  decomposes  along  the  surface  of  the  bubbles  leaves  the  droplet)  when  cp  >  o, . 
Figures  3  and  5  show  that  the  results  of  the  (t,  0)  and  (rj,  £)  models  are  in  good  agreement 
(with  the  exception  of  the  oscillations  in  rs)  in  the  limit  where  0C  — >  1.  The  radial  profiles  of 
temperature  for  the  (r],  £)  model  with  base  case  data  and  0C  =  0.8  are  similar  to  Fig.  4  except 
0max  increases  to  about  7.0  near  the  droplet  surface  when  the  droplet  is  almost  completely 
decomposed.  Although  the  (?/,  £)  model  agrees  well  with  the  (r,  0)  model  for  0C  =  0.99  and 
also  predicts  a  value  of  that  agrees  fairly  well  with  Eq.  (28),  bubbles  probably  connect 
before  0  reaches  0.99.  Therefore,  the  following  results  are  presented  for  the  <;//.  £)  model  with 
oc  =  0.8. 

Figure  6  shows  the  nondimensional  rate  of  mass  loss  (divided  by  A~)  at  the  droplet 
surface,  ips  =  —d0s/dri,  versus  time  for  the  base  case.  As  discussed  above,  the  droplet 
loses  mass  due  to:  gasification  of  the  droplet  surface  (proportional  to  .4]>.  decomposition 
of  bubbles  that  connect  with  the  droplet  surface  (including  the  bubbles  in  the  subsurface 
layer  of  thickness  8  =  Rb),  and  gas  loss  due  to  droplet  surface  regression.  The  gasification 
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Figure  5:  Variation  of  void  fraction  just  below  the  droplet  surface  (/=IDROP)  for  (r.  tr)  and 
(?7,  £)  models  using  data  for  the  base  case.  Different  values  of  or  are  shown  for  the  (//.  £) 
model. 


-  Droplet  Surface 

- Decomposition  in  Connecting  Bubbles 

- Gas  Escape  (surface  regression) 


Figure  6:  Nondimensional  rate  of  mass  loss  (divided  by  4tt)  at  the  droplet  surface.  t.\,  == 
— dips/ dr],  versus  time  for  the  base  case  with  4>c  =  0.8. 
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rate  of  the  droplet  surface  varies  by  less  than  one  order  of  magnitude  throughout  the  droplet 
lifetime.  Droplet  mass  loss  due  to  gasification  from  bubbles  that  connect  with  the  droplet 
surface  exceeds  the  droplet  surface  gasification  rate  at  77  =  7.0,  which  is  approximately  14’ { 
of  the  droplet  lifetime.  The  mass  loss  rate  due  to  bubbles  connecting  with  the  droplet  surface 
increases  sharply  when  o(t=IDROP)  reaches  (pc  (at  t]  =  12.4).  and  exceeds  the  droplet  surface' 
gasification  rate  by  more  than  three  orders  of  magnitude  for  ?/  >  0.5%,*. 

Decreasing  0  [through  Aig  in  Eq.  (20)]  below  unity  increases  bur  does  not  change 
qualitatively  the  radial  profiles  of  9,  R<,.  etc. 

Eqs.  (28)  and  (6)  predict  that  a  p2^.  The  numerical  results  indicate  that 
increases  more  weakly  with  ambient  pressure.  For  p'  =  0.2  atm,  we  find  x  p01-  whereas 
oc  p026  for  p'  —  10  atm.  For  all  the  parameters  studied,  the  temperature  peaks  near  the 
droplet  surface  after  the  droplet  has  decomposed  appreciably.  As  shown  by  Figure  7b.  the 
gas-phase  temperature  decays  more  gradually  with  r  at  lower  p' . 

Effects  of  Other  N ondimensional  Parameters 

The  nondimensional  equations  reveal  that,  the  most  relevant  nondimensional  parameters 
in  this  problem  are  N/m,  k d.  Qd ,  A'/AJ,  and  Ed/{Tl'uT'0).  The  nondimensional  parameter 
TQ/V0  was  not  investigated  for  the  case  of  no  gas-phase  oxidation.  Radial  temperature 
profiles  corresponding  to  the  times  at  which  50%  of  the  initial  liquid  mass  has  decomposed 
and  the  decrease  in  droplet  mass  relative  to  the  initial  droplet  mass  as  a  function  of  time 
are  shown  in  Figures  7  and  8.  respectively,  for  the  above  parameters.  Figure  8  shows  that, 
for  all  cases  except  Qd  —  1  and  kd  =  2.5  x  10-5.  the  droplet  gasifies  very  quickly  after 
approximately  20%  of  the  initial  liquid  mass  has  decomposed. 

The  number  of  bubbles  per  unit  mass,  N/m,  directly  affects  the  void  fraction  [Eq.  (23)]. 
Results  indicate  that  decreases  with  increasing  N/m.  Consistent  with  Eq.  (28).  is 
proportional  to  approximately  ( N/m)~l /3  for  N/m  >  108.  For  much  smaller  N/m  (e.g.. 
N/m  =  104),  r]d,*  oc  (N/m)~023,  indicating  a  small  deviation  from  the  chemical  limit.  As 
N/m  increases,  the  peak  temperature  corresponding  to  the  times  at  which  50%  of  flic  ini¬ 
tial  liquid  mass  has  decomposed  increases  and  the  temperature  distribution  in  the  droplet 
becomes  less  uniform  because  less  time  is  available  for  heat  to  conduct  through  the  droplet 
iFig.  8b). 

The  nondimensional  rate  constant.  kd,  can  be  viewed  as  a  modified  Damkohler  number 
because  it  is  proportional  to  the  ratio  of  the  decomposition  rate  constant  to  the  diffusion  rate. 
Although  decomposition  occurs  in  the  liquid  phase  while  the  diffusivity  in  the  denominator 
of  kd  is  for  the  gas  phase,  varying  kd  is  equivalent  to  varying  the  decomposition  rate  for 
fixed  values  of  the  other  parameters.  Consistent  with  the  scaling  analysis  [Eq.  (28)].  the 
decomposition  rate  is  very  sensitive  to  kd.  For  large  N/m  and  relatively  large  kd  (e.g..  kd  > 
10-4),  the  numerical  model  predicts  that  %,*  is  approximately  proportional  to  kff(N/m)~l  :i 
in  agreement  with  Eq.  (28).  If  //<*.*  oc  kd{N/m)b,  then  defining  «  =  —  (a+3b).  the  dimensional 
droplet  lifetime  varies  with  initial  droplet  radius  as  ifdt  oc  r'2  ~n.  The  chemical  limit  [Eq.  (28)] 
is  a  —  2  while  the  diffusion  limit  is  a  =  0.  Figure  9  summarizes  the  effect  of  k<t  and  N/m 
on  the  dependence  of  droplet  lifetime  ifd  t  on  the  initial  droplet  radius.  The  diffusion  limit 
a  =  0  is  not  predicted  by  the  present  model  because  the  droplet  is  assumed  to  lost'  mass 
solely  due  to  decomposition  and  not  due  to  a  phase-equilibrium  relationship. 
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-  base  case,  n=41.5 

.  kd=2.5xicr\  n=1 19-3 

- -  kd=5.0  x  10”4,  r|=10.3 

- £^^=6.9,11=1 22.9 


Radius,  r 


- base  case,  r|=4l  .5 

.  N/m=105,  r|=75.2 

- -  N/m=107,  ri=21.8 

- p’=10  atm,  r|=74.8 

—  -  •  p’=0.2  atm,  11=34.1 


(a)  (b) 

— .  base  case,  rj=41.5 

. .  Qd=1 ,  r)=91 .4 

- -  Qd=10,  r|=27.3 

- yV=o.o2,  n=4i.o 


(c) 


Figure  7:  Effect  of  parameters  on  the  radial  profiles  of  temperature  (6C  =  0.8)  at  times 
corresponding  to  50%  droplet  mass  loss.  The  base  case  is  kd  =  10-4,  Q(i  =  3.  N/m  =  10'’. 
X'/X[  =  0.2,  p'  =  1  atm,  and  E'd/{H'UT^  =  4.6. 
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1:  base  case 
2:  kd  =  2.5  x  1(T5 
3:  kd  -  ox  1CT5 
4:  E'JCR'Jl)  =  6.9 
5:  E'd/(n’X)  =  2.3 
6:  N/m  =  10s 


1:  base  case 
2:  Qd  =  1 
3:  Q,i  =  10 
4:  //  =  10  atm 
5:  p'  —  0.2  atm 
6:  A9/A ,  =  20 


Time,  ri 


(a) 


(b) 


Figure  8:  Decrease  in  droplet  mass  relative  to  the  initial  droplet  mass.  —Awt,  as  a  function 
of  time  for  the  parameters  k^.  E'd/(TZ'uTg),  X/m,  Qd,  p'.  and  A'/A \. 
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- N/m=108 

. -  Chemical  Limit 


1.0  1 
10 


•5 


10‘ 


Figure  9:  Dependence  of  droplet  lifetime  7]'d  t  on  the  initial  droplet  radius  r'0;  )fll  t  a  r'0‘  " 
The  chemical  limit  is  q  =  2  while  the  diffusion  limit  is  a  =  0. 
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As  ha  decreases,  the  peak  temperature  0max  decreases  and  9  becomes  more  uniform 
in  the  droplet  (Fig.  7a).  Like  the  parameter  kj,  the  nondimensional  activation  energy. 
E'd/(1VUT'0),  directly  affects  the  gasification  rate.  The  droplet  lifetime  depends  more  strongly 
on  E'd/{n'uV0)  for  smaller  values  of  E'd/ (TZ'UT'0)  (i.e.,  closer  to  the  chemical  limit).  For 
<pc  =  0.8,  the  dependence  of  on  E'dj {'R'UT'0)  varies  between  //,/.*  x  and 

Pd,*  oc  e°-oEd^K'‘T’>)  for  2.3  <  E’d/ (1Z'UT'0)  <  6.9.  For  (pc  —  0.99.  this  dependence  varies  between 
Pd,*  oc  eQMEd! ^uTo)  and  pd,t  x  e0'76£”f/(TC“7°).  Although  decreasing  Ed/( H'„T'0)  or  increasing 
kd  increases  rjd t*.  the  temperature  profile  for  a  relatively  small  value  of  E'JCR'^TI,)  =  2.3  is 
relatively  uniform  and  similar  to  the  case  of  a  value  of  kd  that  is  smaller  than  the  base  ease* 
value  (Fig.  7a).  For  smaller  kd,  the  more  uniform  temperature  profile  is  due  to  the  increased 
role  of  conduction  whereas  for  smaller  E'd/(TZ'uTg)  the  more  uniform  9  profile  is  due  to  the1 
lower  activation  temperature  required  for  appreciable  chemical  reactions. 

According  to  Eq.  (28),  Qd  affects  the  characteristic  time  for  droplet  decomposition  only 
indirectly  through  the  evaluation  of  pg  and  9 *.  As  Qd  increases,  gasification  occurs  at  a  higher 
temperature  and  also  pg  decreases.  Substituting  pg  x  Odl  and  a  value  of  9t  greater  than 
unity  into  Eq.  (28)  would  result  in  a  dependence  of  pd.,  on  Qd  stronger  than  //,/.,  x  Qd2 "!. 
Decomposition,  however,  occurs  very  rapidly  once  (p  — ^  o(  near  the  droplet  surface.  Through¬ 
out  the  liquid,  0/  <  1  for  the  majority  of  the  time  before  the  droplet  is  fully  decomposed  (see 
Fig.  4).  Therefore,  the  model  predicts  that  pd.*  has  a  weaker  dependence  on  Qd  (approxi¬ 
mately  p d,*  x  Qd1'2)-  For  all  the  parameters  investigated.  0max  «  Qd  at  the  time  when  ~AYX 
of  the  initial  liquid  mass  has  decomposed  (Fig.  7). 

The  gas-phase  thermal  conductivity,  A',  appears  in  the  parameter  Leg.  which  is  not  varied 
in  the  parametric  study.  Therefore,  changing  the  ratio  of  the  thermal  conductivities.  A'/AJ. 
is  essentially  the  same  as  varying  AJ.  Increasing  A'/A{  (e.g.,  by  decreasing  the  liquid-phase 
thermal  conductivity)  increases  slightly  the  overall  time  for  decomposition  (Fig.  8b)  because 
heat  conducts  through  the  droplet  more  slowly.  Increasing  A'/AJ  also  results  in  a  more  peaked 
temperature  profile  near  the  droplet  surface  (Fig.  7c)  because  of  the  reduced  heat  conduction 
in  the  droplet. 


B.  Results  With  Gas-Phase  Oxidation 

The  ignition  delay  time  and  flame  propagation  rate  are  influenced  by  the  gasification 
rate,  which  is  inversely  proportional  to  ?? d,*.  Hence  it  will  be  shown  that  Eq.  (28).  derived 
above  for  the  chemical  rate  control  limit  of  the  droplet  lifetime  without  gas-phase*  oxidation, 
also  correlates  the  effect  of  decomposition  parameters  on  the  flame  behavior. 

The  flame  location  herein  is  defined  as  the  point  of  maximum  uy./r.  Alternative*  definitions 
such  as  the  point  of  maximum  Yp  do  not  change  the  position  of  the  flame  significantly. 
Unless  otherwise  specified,  all  of  the  calculations  reported  herein  used  kd  =  10 — 1 .  Qd  =  3. 
N/m  =  106.  A'/AJ  =  0.2,  and  E'd/ {VJUT'0)  =  4.6.  Most  of  the  nondimensional  parameters 
were  varied  by  at  least  one  order  of  magnitude  from  their  base  case  values  in  the  parametric 
study  below.  The  following  results  were  recently  submitted  for  archival  publication  (Schiller 
et  al. ,  1997). 
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Flame  Radius,  Flame  Temperature,  and  Ignition  Delay 

In  classical  droplet  burning,  the  ratio  Rf\/Rs  asymptotically  tends  to  a  constant,  and  the 
flame  radius  first  increases  but  then  decreases  due  to  the  shrinking  of  the  droplet.  Similar 
behavior  is  predicted  by  the  present  model  if  decomposition  is  not  considered  and  classical 
droplet  vaporization  boundary  conditions  are  used  at  the  droplet  surface.  Conditions  at  the 
interface  in  the  model  with  decomposition,  however,  differs  from  classical  droplet  models  in 
two  ways:  (1)  the  gasification  rate  is  determined  from  the  decomposition  model  rather  than 
from  the  use  of  a  Clausius-Clapeyron  relation  combined  with  a  Stefan  flow  assumption,  and 
(2)  the  gasification  process  is  exothermic.  As  a  result,  the  model  predicts  that  the  variation 
of  flame  radius  with  time  is  quite  different  than  that  of  classical  droplet  combustion. 


-  Base  Case  (p=1  atm,  N/m=106) 

.  p=10atm 

- N/m=103 


:  .  i 

300  400 


Figure  10:  Flame  radius  and  flame  temperature  versus  time  for  different  values  of  ambient 
pressure  and  number  of  bubbles  per  unit  mass.  The  droplet  radius  (multiplied  by  10)  is  al>o 
shown  for  the  base  case.  The  droplet  lifetime,  is:  63.4  (base  case),  59.5  (;/  =  10  atm*, 
and  333.7  ( N/m  =  103). 


Over  the  investigated  range  of  the  fundamental  parameters,  R{\  first  increases  nearly 
linearly  with  r/  after  an  ignition  delay  period  (Fig.  10).  The  ignition  delay  time  increases  with 
decreasing  gasification  rate  (e.g..  due  to  a  decrease  in  A^  or  N/m),  in  a  manner  qualitatively 
similar  to  Eq.  (28).  Plots  of  Ra  and  6a  versus  ??  for  k <<  =  10"°  and  N/m  =  10°  are  nearly 
identical  to  those  shown  in  Fig.  2  for  Av  =  10-4  and  N/m  =  103.  Eq.  (28)  predicts  that  a  one 
order  of  magnitude  decrease  in  A’ a  —  has  the  same  effect  on  the  droplet  lifetime  as  a  three 
order  of  magnitude  decrease  in  N/m.  The  flame  behavior  thus  has  a  qualitatively  similar 
dependence  on  the  parameters  that  affect  the  decomposition  rate. 

As  i]  increases,  the  droplet  radius  first  increases  due  to  bubble  growth,  but  then  su!>- 
sequently  decreases  as  the  droplet  gasifies  (Fig.  10).  R&/Rs  also  increases  with  time*  for 
i]  <  Vd,*-  Unlike  classical  droplet  burning,  no  simple  relationship  exists  for  the  variation  of 
with  time.  After  decomposition  is  artificially  terminated  at  //  =  ly...  R\\  increases 
but  then  subsequently  decreases  with  //.  similar  to  classical  droplet  combustion.  Tin'  flame 
radius  might  vary  differently  with  time  while  the  droplet  decomposes  as  compared  to  classi¬ 
cal  droplet  combustion  because  the  gasification  rate  due  to  decomposition  is  more  than  one 
order  of  magnitude  greater  than  that  due  to  classical  vaporization.  Apparently,  the  higher 
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gasification  rate  causes  the  flame  to  move  away  from  the  droplet  surface  more  rapidly. 


- Qc=20 

.  Qc=58  (Base  Case) 

- Qc=100 

—  -  No  oxidation 


Figure  11:  Flame  temperature  and  droplet  surface  temperature  versus  time  for  different 
values  of  Qc.  The  droplet  surface  temperature  is  also  shown  for  the  case  of  no  gas-phase 
oxidation  ( kc  =  Qc  =  0).  %.*  varies  between  63  and  64  for  these  cases. 


Over  the  investigated  range,  parameters  that  directly  affect  only  the  decomposition  rate 
(e.g.,  N/m.  kd,  etc.)  or  the  rate  of  oxidation  do  not  appreciably  change  the  maximum 
temperature  in  the  gas  phase.  Figures  10  and  11  show  that  the  flame  temperature  does 
not  vary  significantly  with  time  soon  after  the  flame  is  initiated.  As  expected.  dt \  is  fairly 
insensitive  to  ambient  pressure  and  increases  as  the  amount  of  heat  released  due  to  gas- 
phase  oxidation  is  increased.  For  benzene,  the  adiabatic  flame  temperature  with  a  constant 
value  of  specific  heat  and  Vb2,oc  =  0.23  is  da  =  1  +  0.07(2,.)  which  is  consistent  with  the 
model’s  predicted  values  of  da  for  Qd  —  0  (no  heat  release  due  to  decomposition).  As  shown 
by  Figs.  10  and  11.  the  base  case  value  of  Qd  —  3  causes  da  to  increase  slightly  above 
the  adiabatic  flame  temperature  (e.g.,  by  approximately  6%  for  Qc  =  58).  Also  consistent 
with  the  theoretical  values  of  adiabatic  flame  temperature  for  constant  Cp.  On  increase's  from 
approximately  2.0  to  9.4  when  l’o,.oc  is  changed  from  0.05  to  0.46  (Fig.  12). 
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Figure  12:  Flame  temperature  versus  time  for  three  values  of  1  o,,oc  (curve  labels). 
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Increasing  Qd  from  0  to  5  decreases  the  ignition  delay  time  from  28  to  23  (Fig.  13) 
and  increases  6f\  from  5.2  to  5.7.  As  shown  by  Fig.  13.  higher  values  of  Q,i  lead  to  greater 
gasification  rates  and  thus  faster  flame  propagation.  For  Q,i  =  0,  the  droplet  lifetime  is 
almost  equal  to  the  time  required  for  complete  oxidation  (//,/.*  -  304.8.  //,..*  =  312.4).  whereas 
oxidation  occurs  relatively  long  after  the  droplet  has  fully  decomposed  for  Q,i  =  5  (//,/  .  = 
48.0.  tic,*  =  149.2). 
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Figure  13:  Flame  radius  and  nondimensional  gasification  rate  as  functions  of  time  for  three 
values  of  Qd. 


Effect  of  Gas-Phase  Oxidation  on  Droplet  Lifetime 

Gas-phase  oxidation  does  not  appreciably  affect  the  droplet  lifetime  for:  (1)  the  base  case 
values  of  the  nondimensional  parameters,  or  (2)  any  variation  of  parameters  that  results  in  a 
fast  rate  of  gasification  compared  to  oxidation.  The  insensitivity  of  ?/(/.*  to  gas-phase  oxidation 
occurs  because  the  droplet  lifetime  is  governed  primarily  by  the  chemical  rate'  control  limit 
of  liquid-phase  decomposition.  The  temperatures  in  the  droplet  rise  quickly  due  to  the  heat 
released  by  decomposition  before  they  are  affected  by  heating  due  to  conduction  from  t In¬ 
flame.  Gas-phase  oxidation  does,  however,  appreciably  affect  tpt,*  if  the  rate  of  decomposition 
is  slowed  significantly,  thus  giving  more  time  for  conduction  from  the  flame  to  affect  tin- 
droplet  heating.  Figure  14  shows  the  nondimensional  total  mass  rate  of  gasification  and 
oxidation  as  functions  of  time  for  two  values  of  N/m.  For  N/m  =  10()  (base  case  value), 
the  droplet  is  gasified  at  =  63.4.  whereas  oxidation  is  not  complete  until  //<.,,  =  1GC.2. 
Because  gasification  is  relatively  fast  compared  to  oxidation  for  the  base  case.  gas-phase 
oxidation  decreases  r id,*  bv  less  than  19c.  Conversely,  for  X/m  =  10'}.  decomposition  occurs 
more  slowly,  thus  giving  more  time  for  conduction  from  the  flame  to  affect  the  droplet  heating. 
Unlike  the  base  case,  at  any  given  time  the  total  mass  that  has  oxidized  is  comparable  to 
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that  which  has  gasified  for  N/m  =  103.  In  this  case,  gas-phase  oxidation  decrease's  //,/.,  by 
16%.  Plots  of  gasification  and  oxidation  rates,  as  well  as  the  effect  of  gas-phase  oxidation  on 
rid,*-  for  kd  =  10"5  and  N/m  —  103  are  similar  to  those  for  kd  =  10~4  and  N/m  =  1()3. 


n 

Figure  14:  Total  mass  rate  of  gasification  and  oxidation  as  functions  of  time  for  two  value's 
of  N/m. 

As  the  ambient  pressure  increases,  both  the  droplet  lifetime  and  the  time  required  for 
complete  oxidation  decrease  (Fig.  10),  and  the  ratio  i)c.*/>}d.*  approaches  unity.  In  general.  the' 
effect  of  gas-phase  oxidation  on  the  droplet  lifetime  decreases  as  *//;</.*  becomes  greater 
than  unity.  For  example,  the  decrease  in  //<*.*  due  to  gas-phase  oxidation  is  19%  for  //  = 
20  atm  (ric^/rjd, *  =  1.09),  7.5%  for  //  =  10  atm  {i]c,*/rid.*  =  1.23).  and  0.6%  for  //  =  1  atm 
(Vc,*/Vd.*  =  2.62). 

The  droplet  lifetime  is  insensitive  to  large  changes  in  Ay.  E'c/{7l'uT'0).  QCl  and  Yq.2,0 c.  For 
example,  i]d ,»  decreases  by  roughly  1%  when  either  Qc  is  increased  from  20  to  100,  Ya.,^  is 
increased  from  0.05  to  0.46,  or  kc  is  increased  by  two  orders  of  magnitude.  The  insensitivity  to 
changes  in  the  oxidation  rate  constants  is  indicative  of  a  fast  reaction.  Although  Qv  and  1  A,.x 
affect  the  flame  temperature,  they  do  not  appreciably  change  for  base  case  values  of  the 
other  parameters  because  the  droplet  lifetime  is  governed  by  liquid-phase  decomposition.  Tin1 
plot  of  droplet  surface  temperature  (9S)  versus  time  in  Fig.  11  indicates  that  the  increase  in 
flame  temperature  with  Qc  does  increase  6S  for  much  of  the  droplet  lifetime.  But  the  droplet 
temperature  profile  and  thus  the  gasification  rate  are  affected  more  by  the  exothermic  rate 
of  heat  release  due  to  liquid  decomposition  than  by  the  increased  rate  of  heat  conduction 
from  the  flame.  The  insensitivity  of  6S  (and  0  throughout  the  droplet)  to  6 a  (through  Qr)  is 
especially  evident  in  the  last  15%  of  the  droplet  lifetime,  when  liquid-phase  decomposition 
dominates  the  droplet  heating. 

As  shown  by  Fig.  15,  the  droplet  almost  fully  decomposes  before  one-third  of  the  fuel 
has  oxidized.  Higher  ambient  temperature  results  in  more  rapid  heat  conduction  to  the 
droplet  and  more  rapid  exothermic  gasification  (decomposition)  with  subsequent  further 
liquid  heating.  The  droplet  lifetime  therefore  decreases  with  an  increase  in  T'^/T'0.  [T'0  was 
fixed  at  300  K  for  these  simulations.)  Gas-phase  oxidation,  however,  does  not  affect  //(/i 
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Figure  15:  Total  mass  of  fuel  ( rnF )  oxidized  and  decomposed  as  functions  of  time  for  four 
values  of  T'^/T'g. 

over  the  investigated  range  of  ambient  temperatures  (1.1  <  T'^./T'0  <  C.67).  Moreover.  the 
rapid  increase  of  the  flame  radius  (Fig.  10)  increases  the  time  required  for  heat  to  conduct 
from  the  flame  to  the  droplet  surface.  Hence  the  droplet  surface  temperature  increases  more 
rapidly  with  increasing  T^/T^  than  with  increasing  flame  temperature  (compare  Fig.  10 
with  Fig.  11).  The  heat  release  due  to  liquid-phase  decomposition  causes  0,  to  exceed  the 
flame  temperature  (Figs.  11  and  16).  As  shown  by  Fig.  16.  the  droplet  surface  temperature 
decreases  rapidly  with  time  soon  after  the  liquid  stops  decomposing.  For  T'^/T^  =  1.5.  the 
flame  extinguishes  before  30%  of  the  initial  fuel  mass  is  oxidized.  Ignition  does  not  occur 
for  VJT0  <  1.2. 


Figure  16:  Flame  temperature  and  droplet  surface  temperature  as  functions  of  time  for 
T'x/Tg  =  1.5  (%,,  =  180.6,  flame  extinguishes)  and  T'^jT'0  =  2.0  (//,/,*  =  106.0.  //,..*  -  174.4). 
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Gas-Phase  Profiles 

As  stated  above,  the  droplet  lifetime  is  affected  by  gas-phase  oxidation  if  the  rate  of  liquid- 
phase  decomposition  is  decreased  sufficiently.  The  droplet  lifetime  is  reduced  by  27 due 
to  gas-phase  oxidation  for  N/m  =  102.  As  shown  by  Fig.  17a.  there  is  appreciable  heating 
of  the  droplet  by  gas-phase  conduction  from  the  flame  over  much  of  the  droplet  lifetime 
because  %,*  is  relatively  large.  That  is,  the  time  scale  for  heat  conduction  from  the  flame  to 
the  droplet  becomes  comparable  or  even  smaller  than  that  for  liquid-phase  decomposition. 
Comparatively,  for  N/m  =  106,  the  heat  from  the  flame  does  not  have  time  to  conduct  to 
the  droplet  surface  before  77  =  As  shown  in  Fig.  17b.  heat  actually  conducts  away  from 
the  droplet  surface  during  most  of  the  droplet  lifetime  for  N/m  =  10°.  Results  indicate  that 
the  speed  of  the  flame.  clRa/di],  decreases  with  rid.*-  The  time  required  for  heat  to  conduct 
from  the  flame  to  the  droplet  surface  varies  with  y/Rn  —  1.  Therefore,  the  reduction  in  flame 
speed  as  r/d.*  is  reduced  contributes  to  the  decrease  in  the  time  scale  for  heat  conduction  from 
the  flame  to  the  droplet.  Similar  to  classical  droplet  burning,  the  profiles  of  Y,  in  Fig.  17 
indicate  that  the  flame  is  diffusion  controlled,  with  almost  no  oxygen  in  the  region  r  <  7?f  1 
and  no  fuel  vapor  in  the  region  r  >  Rr. 

CONCLUSIONS 

A  simplified,  spherically-symmetric  model  of  the  transient  hearing  of  an  energetic  liquid 
fuel  droplet  has  been  developed.  The  decomposition  of  the  energetic  fuel  was  examined 
for  a  range  of  nondimensional  parameters  including  a  modified  Damkdhler  number,  heat  of 
decomposition,  activation  energy,  number  of  bubbles  per  unit  mass,  and  ratio  of  thermal 
conductivities. 

Without  the  effects  of  gas-phase  oxidation,  detailed  results  of  the  base  case,  which  were 
qualitatively  similar  to  all  of  the  simulations/revealed  the  following.  The  bubble  radius  grows 
with  radial  position  and  time,  causing  the  droplet  to  swell  initially.  A  model  that  incorporates 
decomposition  of  the  droplet  surface  predicts  that  the  droplet  radius  continuously  increases 
with  time  until  surface  decomposition  becomes  appreciable,  after  which  time  tin*  droplet 
mass  and  droplet  radius  decrease  continuously  with  time.  The  time  required  for  the  droplet 
to  decompose  decreases  as  4>c  is  increased  or  ft  is  decreased  in  the  model  that  includes  droplet 
surface  decomposition.  A  model  that  assumes  the  droplet  mass  decreases  due  to  bubbles 
discontinuously  bursting  at  the  droplet  surface  predicts  that,  after  0  — »  1  just  below  the 
droplet  surface,  the  intermittent  bursting  of  bubbles  into  the  gas  film  causes  small  oscillations 
of  the  droplet  radius  as  it  intermittently  swells  (due  to  sensible  heating  of  the  bubbles)  and 
shrinks  (due  to  the  bubble  bursting).  The  mean  value  of  r,  decreases  with  the  mass  of  the 
droplet.  The  results  of  the  two  models  are  in  good  agreement  (with  the  exception  of  the 
oscillations  in  rs)  in  the  limit  where  0C  — »  1.  The  temperature  peaks  just  outside  of  t Ik* 
droplet  surface  due  to  decomposition.  Heat  conduction  is  dominant  during  the  early  stages 
of  droplet  heating,  then  the  chemical  energy  release  dominates  later. 

Consistent  with  simplified  scaling  for  the  limit  of  chemical  rate  control,  the  characteristic 
time  for  decomposition  (rjd,*)  is  strongly  dependent  011  kd  and  E'd/{TVaT'0).  and  less  strongly 
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Figure  17:  Gas-phase  profiles  of  9  and  }  }  at  different  values  of  i]  and  .V/ in 


28 


on  iV/m,  p',  and  Qd.  Increasing  A'/AJ  increases  slightly  and  results  in  a  more  peaked 
temperature  profile.  Over  the  range  of  parameters  investigated.  ifd  „  is  proportional  to  r'„ 
to  a  power  between  0.1  and  0.5,  indicating  that  the  chemical  kinetic  control  of  the  droplet 
lifetime  is  stronger  than  the  diffusion  control. 

With  the  effects  of  gas-phase  oxidation,  the  flame  radius  is  predicted  to  increase  nearly 
linearly  with  time  before  subsequently  decreasing.  This  behavior  differs  from  classical  droplet 
burning  due  to  the  exothermic  decomposition  process  that  determines  the  gasification  rate. 
Because  the  ignition  delay  time  and  flame  propagation  rate  depend  on  the  gasification  rate 
of  the  droplet  (typically  controlled  by  liquid-phase  decomposition  rather  than  gas-phase* 
conduction),  simplified  scaling  previously  derived  for  the  droplet  lifetime  also  correlate's  the* 
effect  of  decomposition  parameters  on  the  flame  behavior.  Decreasing  the  decomposition 
rate  (e.g.,  by  reducing  kd  or  N/m )  increases  the  ignition  delay  time  and  reduces  the*  flame' 
propagation  rate. 

For  the  selected  base  case,  liquid-phase  decomposition  increases  the  flame  temperature 
by  approximately  6 Si.  Therefore,  energetic  liquid  fuels  can  be  expected  to  vaporize  and 
burn  more  rapidly  and  to  yield  more  thermal  energy  than  conventional  hydrocarbon  fuels. 
The  droplet  lifetime  for  the  base  case,  however,  is  not  affected  significantly  by  gas-phase 
oxidation  because  the  droplet  fully  decomposes  before  heat  conduction  from  the  flame  can 
appreciably  affect  the  droplet  heating.  As  the  decomposition  rate  is  reduced,  the  time  scale 
for  heat  conduction  from  the  flame  to  the  droplet  becomes  comparable  to  that  of  liquid 
decomposition,  and  hence  gas-phase  oxidation  significantly  reduces  In  general,  the 

effect  of  gas-phase  oxidation  on  the  droplet  lifetime  decreases  as  becomes  greater 

than  unity. 

Gas-phase  profiles  of  0  and  I,  for  the  base  case  (A/m  =  106)  indicate  that:  (1)  the  flame 
is  diffusion  controlled.  (2)  the  flame  radius  increases  with  /;.  and  (3)  heat  from  the  flame 
does  not  conduct  into  the  droplet  over  most  of  the  droplet  lifetime.  If.  however,  A  / m  is 
reduced  to  102,  thus  slowing  the  decomposition  rate,  then  there  is  sufficient  time  for  heat  to 
conduct  from  the  flame  to  the  droplet.  In  this  case,  gas-phase  oxidation  reduces  the  droplet 
lifetime  by  approximately  27%. 
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